Magnetic interference patterns in O-vr SIFS Josephson junctions: effects of asymmetry 

between and tt regions 
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We present a detailed analysis of the dependence of the critical current Ic on the magnetic field B 
of 0, TT, and O-tt superconductor- insulator- ferromagnet-superconductor Josephson junctions. Ic{B) 
of the and tt junction closely follows a Fraunhofer pattern, indicating a homogeneous critical 
current density jdx). The maximum of Ic{B) is slightly shifted along the field axis, pointing to a 
small remanent in-plane magnetization of the F-layer along the field axis. Ic{B) of the O-vr junction 
exhibits the characteristic central minimum. 7c however has a finite value here, due to an asymmetry 
of jc in the and tt part. In addition, this Ic{B) exhibits asymmetric maxima and bumped minima. 
To explain these features in detail, flux penetration being different in the part and the tt part 
needs to be taken into account. We discuss this asymmetry in relation to the magnetic properties 
of the F-layer and the fabrication technique used to produce the O-tt junctions. 

PACS numbers: 74.50.-|-r,85.25.Cp 74.78. Fk 74.81.-g 

Keywords: 0-pi Josephson junction; ferromagnetic Josephson junction; SIFS Josephson junction; fractional 
magnetic flux quantum; semifluxon 



I. INTRODUCTION 

While predicted more than 30 years ago'^'^, due to 
the severe technological requirements, the experimental 
study of TT Josephson junctions became an intense field 
of research only recently. Superconductor-ferromagnet- 
superconductor (SFS) Josephson junctions were success- 
fully fabricated and studied''"''''''''. SFS junctions however 
typically exhibit only very small (metallic) resistances R, 
making this type of junctions less suitable for the study 
of dynamic junction properties as well as for applications, 
where active Josephson junctions arc required. To over- 
come this problem, an additional insulating (I) layer can 
be used to increase R, although at the expense of a highly 
reduced critical current density jc''^''''^'\ 

In a SFS or SIFS junction the proximity effect in the 
ferromagnetic layer leads to a damped oscillation of the 
superconducting order parameter in the F-layer. Thus, 
depending on the thickness dp of the F-layer, the sign 
of the order parameters in the superconducting elec- 
trodes may be equal or not. While in the first case a 
conventional Josephson junction (a "0 junction") with 
Is = Ic sin(/x) is realized, in the latter case a "tt junction" 
is formed where the Josephson current I^ obeys the rela- 
tion Is = Ic sin(/x + tt) = —IcShi(fi). Here Jc > is the 
junction critical current and /x is the phase difference of 
the order parameters in the two electrodes. 



The combination of a and a tt part within a sin- 
gle Josephson junction leads to a "O-tt" Josephson junc- 
tion. Depending on several parameters of the and 
TT part, a spontaneous fractional vortex may appear at 
the O-TT boundary In case of long junctions with 
length L 3> Aj the vortex contains a flux equal to a 
half of a flux quantum $0 « 2.07 x lO^^^Tm^. Here 
Aj « •\/$o/(47r/iojcAL) is the Josephson length; /io is 
the magnetic permeability of the vacuum and A^ is the 
London penetration depth of both electrodes. 

Up to now three different types of O-tt Joseph- 
son junctions exist. One approach makes use of the 
d^2_y2 wave order parameter symmetry in cuprate 
superconductors'"'^''^ ''^'''^'^'. Another approach is to 
use standard Nb/Al-A^Os/Nb Josephson junctions 
equipped with current injectors'''''^ which allow to cre- 
ate any phase shift. O-tt Josephson junctions were also 
produced (accidentally) by SFS technology^"'-". The 
first intentionally made O-tt SIFS junction including ref- 
erence and TT Josephson junctions fabricated in the 
same run were recently realized^'. Some static and dy- 
namic properties of this type of O-tt junction were stud- 
ied experimentally-''--'- '. Relevant theoretical work on 
SIFS junctions can be found in'"''^^. 

The aim of the present paper is to provide a careful 
analysis of the magnetic field dependence of the junction 
critical current Ic{B) in order to characterize these novel 
type of junctions as accurately as possible. The (short) 
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FIG. 1; (Color online) (a) Sketch of a O-tt SIFS Josephson 
junction with step-like F-layer to create a 0-coupled part (F- 
layer thickness di) and a n part (d2). LqjLtt denote the 
length in each part, (b) Schematic jcidp) dependence for 
SIFS Josephson junctions. For the ferromagnetic thicknesses 
di and d2 = di + AdF the critical current densities jc{di) and 
jc{d2) have similar absolute values {jc{di) — jc{d2) + Ajc). 



junction we discuss has a length L m Xj. As we will see, 
the measured Ic{B) can be reproduced very well when, 
apart from asymmetries of the critical current densities 
in the and tt parts, asymmetric flux penetration into 
the and tt parts is taken into account. 

The paper is organized as follows: In section II the 
SIFS junctions are characterized in terms of geometry, 
and the properties of the F-layer are further characterized 
by measuring the magnetization of a bare Nio.6Cuo.4 thin 
films with thickness comparable to the F-layer used for 
the junctions. In the central section III the magnetic field 
dependence of the critical current of the SIFS junctions 
is discussed. Section IV contains the conclusion. 



II. SAMPLE CHARACTERIZATION 

Fig. 1(a) shows a sketch of the O-tt junction used in 
the experiment. The superconducting bottom and top 
layers consist of Nb with the thicknesses ti = 120 nm 
and t2 = 400 nm, respectively. As for standard Nb tun- 




100 nm 



FIG. 2: (Color online) Optical image (top view) of a 330 x 
30 /im^ window junction with junction and idle regions. 



ncl junctions an AI2O3 layer was used as tunnel barrier. 
Its thickness is dj « 0.9 nm, determined from dynamic 
measurements. For the ferromagnetic layer we use the di- 
luted ferromagnet Nio.6Cuo.4. To form a O-tt junction the 
junction is divided into two parts differing by the thick- 
ness of the F-layer. While in one half of the junction the 
thickness di is chosen such that coupling is realized, in 
the other half the F layer thickness d2 is used to realize 
TT coupling. In order to have approximately symmetric 
junctions, di and d^ should be such that the critical cur- 
rent densities of the two halves are about the same and 
as large as possible, see Fig. 1(b). 

Details of the fabrication technique can be found in 
Refs. 22,20. The main feature is a gradient in the fer- 
romagnetic Nio.6Cuo.4 layer along the y direction of the 
4" wafer, in order to allow for a variety of and tt cou- 
pled junctions differing in their critical current densities. 
In addition, by optical lithography and controlled etch- 
ing, parts of the F-layer are thinned by ls.dp « SA, such 
that coupling is achieved in these parts. Thus the chip 
contains un-etched parts with F-layer thickness dpiu), 
as well as uniformly etched parts with F-layer thickness 
— Adp. Thus, at a fixed y-position we have two dif- 
ferent ferromagnetic thicknesses allowing for patterning 
a set of three junctions: 

• a junction with F-layer thickness di and critical cur- 
rent density j° = jc(di) 

• a TT junction with F-layer thickness ^2 and critical cur- 
rent density = jc{d2) 

• a stepped O-tt junction with thicknesses c?i, d2 and crit- 
ical densities jj?, in and n halves. 

For the values di ~ 5.05 nm and d2 = 5.33 nm we 
achieved j° « 2.1 A/cm^ and w 1.7 A/cni^ at 
T = 4.2 K, as estimated from and tt reference junctions. 
Due to the different temperature dependence of j° and 
(see Ref. 8) these values change to j° ~ ~ 2.2 A/cm^ 
at T = 2.65K. 

All junctions had the same geometrical dimensions 
330 X 30 firn^, see Fig. 2. The superconducting electrodes 



3 




-400 -300 -200 -100 100 200 300 



B(mT) 

FIG. 3: (Color online) Magnetization curves of a Nio.6Cuo.4 
thin film with 10 nm thickness at T = 5 K probed by SQUID 
magnetometry. The magnetic field was applied either in-plane 
(squares) and out-of-plane (circles). The inset shows a mag- 
nification at small magnetic fields. 



extend well beyond the junction area, leading to an idle 
region around the junction affecting the Josephson length 
Aj. Ignoring this correction, using — — 2.2 A/cm^, 
as measured at T = 2.65 K, one finds Aj « 260 /im, 
i.e. L w 1.2Aj as in Ref. 21. The idle region of width 
Wi^i + Wi^2 in y direction leads to an effective Josephson 
length- ' 

A,,,ff = A,yH- -, 

with the junction width and the inductances (per 
square) of the superconducting films forming the junc- 
tion electrodes fiod'j and the idle regions fiod'j. For 
our junction we get Xj^cS — 1-7A,/, with W — 30 /im. 
Wis + Wi,2 = 100nm,'d'j = 194 nm, and d'j ^ 350 nm. 
Therefore the normalized junction length at T = 2.65 K 
is / = L/Xj^cS ~ 0.76 and we clearly are in the short 
junction limit. 



Magnetic properties of the F-layer 

In order to investigate the magnetic properties of 
the Nio.6Cuo.4 alloy used for the F-layer we performed 
measurements of the magnetization via SQUID magne- 
tometry. The sample was a 10 nm thin Nio.6Cuo.4 film 
deposited directly on a Si02 substrate. The obtained 
magnetization curves (after diamagnetic correction) at 
T = 5K are shown in Fig. 3 for the magnetic field 
applied in-plane or out-of plane. The magnetic moments 
for the out-of plane and in-plane component saturate 
at almost equal m « 6.5 x 10~^ emu corresponding to 
a saturation magnetization M = 130 kA/m. Using the 
density p = 8.92 g/cm^ (bulk value) of the F-layer and 



the molar weight 60.6 g we can estimate the atomic 
saturation moment mat = 0.16/i5, in good agreement 
with TTiat = 0.15 fiB found in literature"'^. 
In the inset of Fig. 3 the hysteresis of the magnetiza- 
tion curves is shown at small applied magnetic fields. 
Rcmanence can be seen for the in-plane as well as the 
out-of plane curves. The inversion of the magnetizations 
is smooth, indicating a multiple domain state. The 
magnetic field necessary to fully magnetize the magnetic 
film in-plane is in the order of 10 mT, whereas the out-of 
plane magnetization saturates above about 100 mT. 
Therefore we expect the in-plane magnetization to be 
energetically favorable. 

Both saturation fields are orders of magnitude larger 
than the in-plane fields typically used for SIFS critical 
current versus magnetic field measurements. In the 
following we estimate an upper limit by how much 
the Ic{B) pattern (of a junction or a tt jimction) 
would shift along the field axis for an in-plane, fully 
saturated ferromagnetic layer. Our measured saturation 
magnetization Al = 130 kA/m yields a magnetic 
induction ^qM ~ 0.163 T. A cross section of length L 
and a thickness dp encloses an intrinsic magnetic flux 
= dp ■ L ■ p^qM. For L = 330 /im and dp = b nm 
the magnetic flux is $a/ = 129 x <i>o. Thus, the Ic{B) 
pattern would be shifted along the field axis by about 
129 periods, while in experiment typically shifts of 
much less than one period are observed. Further, nearly 
all our SIFS junctions had mirror-symmetrical Ic{B) 
patterns for |_B| < 1 mT, again strongly indicating that 
the F-laycr is in a multiple domain state with a very 
small in-plane net magnetic flux"'\ The out-of-plane 
net magnetic flux has to be small too. As we will see 
in the next section, for the and tt junctions highly 
symmetric Ic{B) patterns can be measured. If the 
out-of-plane magnetic flux were very large, one would 
expect a large number of Abrikosov vortices penetrating 
the superconducting layers, making the Ic{B) of SIFS 
junctions with a planar F-layer strongly asymmetric. 
The ferromagnetic properties of a comparable ferro- 
magnetic compound, Cuo.47Nio.53, were investigated 
recently via anomalous Hall voltage measurements and 
Bitter decoration techniques of the magnetic domain 
structures'*", indicating a magnetic anisotropy and a 
magnetic structure with domains of about 100 nm in 
size. Both Hall and Bitter decoration measurements 
are only sensitive to out-of-plane components of the 
magnetic fields, and the growth conditions of the CuNi 
sample in Ref. may infiuence its magnetic properties. 
Nevertheless it supports our experimental findings of a 
very small in-plane magnetization for zero field cooled 
samples and a multiple domain state in the F-layer of 
our SIFS devices. 
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III. CRITICAL CURRENT VS. MAGNETIC 
FIELD 

In order to measure the magnetic field dependence of 
the critical currents of our junctions, the samples were 
mounted in a glass-fiber Helium cryostat surrounded by 
a triple mu-metal shield. To minimize external noise the 
whole setup was placed in a high-frequency screnning 
room, the current leads were low-pass filtered, and all 
electronics within the screnning room was powered by 
batteries. The sample was initially cooled from room- 
temperature down to 4.2 K with the sample mounted 
inside the magnetic shield. To remove magnetic flux 
sometimes trapped in the superconducting electrodes 
the sample was thermally cycled to above the super- 
conducting transition temperature Tc. To determine 
Ic we used a voltage criterion of = 0.5 /.(V. The 
current- voltage {IV) characteristics and Ic{B) were mea- 
sured for all three junctions at various temperatures 
T = 4.2... 2.65 K. The magnetic field B was applied 
along the y direction see Fig. 1(a). 
Figure 4 shows measurements of Ic{B) at (a) T « 2.65 K 
and (b) T = 4.2 K. Together with the experimental data 
we plot theoretical curves using the analytic expressions 
valid for short junctions having homogenous critical cur- 
rent density. 

For the and tt junctions one has the Fraunhofer pat- 
tern: 



/r(i3)==/°-(0) 



sin(^) 



*0 



(1) 



where $/$o = BLA/^q is the number of the applied 
flux quanta through the normalized junction area LA, 
with A = di + dp + Xl tanh(ii/2AL) + Xl tanh(t2/2AL). 
For a symmetric, short O-tt junction the analytical ex- 
pression is given by *^"'^": 



{ <S>\ 



(2) 



At T = 2.65 K the reference junctions have basically 
the same maximum critical current of w 220 fiA and 
« 217 //A and are fitted very well by the standard 
Fraunhofer curve given by Eq. (1). Note that the max- 
imum is shifted along the B axis by a few percent of 
one flux quantum. For reference we also show by a dot- 
ted horizontal line the /c-detection limit /c,min = Vcr/R 
set by the finite voltage criterion. Here R denotes the 
(subgap) junction resistance at small voltage. R was es- 
timated from the corresponding IV-cuives shown in the 
insets of Fig. 4. For the measurements at T = 2.65 K 
this line is marginally shifted from zero. 

Looking at I^~'^{B) of the stepped O-tt junction at 
T = 2.65 K (see bottom graph of Fig. 4(a)), we see 
that the agreement between the analytical expression 
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FIG. 4: (Color online) Ic{B) measurements and theoretical 
curves (short junction model) for 0, tt and O-tt Josephson junc- 
tions at (a) T r; 2.65 K and (b) T 4.2 K. In the top graphs 
of (a) and (b) data for the and tt junctions are shown by 
solid symbols; the Fraunhofer curves Eq. (1) are shown by the 
solid lines. In the lower graphs data for the O-tt junction are 
shown by symbols; the solid line corresponds to Eq. (2). For 
the theoretical curve shown by the dashed lines an asymmetry 
in the critical current densities S = (jj? — jc)/ijc + jc) ~ 0.18 
in (a) and 0.33 in (b) has been assumed. The horizontal dot- 
ted lines show the value of /c,min. The insets in (a) and (b) 
show JV^-curves for all three Josephson junctions, using the 
same symbols as for the Ic{B) patterns. 



Eq. (2) and the measurement is worse than for the ref- 
erence junctions. For example the central minimum of 
I^~'^{B) is reproduced qualitatively, however, apart from 
a slight shift to positive magnetic field values, it does 
not reach zero critical current and is U-shaped in con- 
trast to the V-shaped central minimum predicted by 
Eq. (2). Further, the side maxima in I^~'^{B) at the 
magnetic field ztBm = 2^q/LA are below the theoretical 
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value of 0.72/^. Additionally we found a small asymme- 
try of the maxima of 4 %, i.e. « 0.66 
and I^"" {+Bm) / 1° « 0.64. Finally, the first side min- 
ima of I^~'^{B) were reached at the same magnetic field 
(<f>/$o = ±2) as the second minima of the Ic{B) of the 
reference junctions, but exhibit bumps and do not reach 
zero-level defined by the /c,min hue. 

All discrepancies to the calculated pattern, especially 
the non-vanishing minima, are not due to our measure- 
ment technique. All characteristic features are well above 
our /c detection limit, drawn by the dotted line in the 
bottom graph of Fig. 4(a). The U-shaped central mini- 
mum /c''^(0) could be due to fluctuations in the apphed 
magnetic field. However, careful measurements using su- 
perconducting magnetic field coils in persistent mode to 
exclude any magnetic field noise showed no further de- 
crease of the minimum. An improved fit can be achieved 
by assuming that the critical current densities of the two 
halves of the O-tt junction arc not identical, i.e. are dif- 
ferent from the respective j° and of the reference 
junctions (e.g. caused by some gradient of the ferro- 
magnetic thickness along x direction; the distance be- 
tween reference and stepped junctions on the chip is 
about 2 mm). The dashed line in the bottom graph 
of Fig. 4(a) shows the result of a corresponding calcu- 
lation (the procedure is discussed further below) using 
6 = {jc-jc)/{jc +jc) = 0-18. While the critical current 
value of the central minimum is reproduced reasonably 
well, the other discrepancies remain. 
Fig. 4(b) shows data for T — 4.2 K. The critical currents 
of the and tt reference junction differ by » 22%, but still 
are reasonably well described by the Fraunhofer pattern 
Eq. (2). The main discrepancy between fit and measure- 
ments can be found in the minima of Ic{B). The exper- 
imental minima do not reach zero current, which at this 
temperature is due to the finite voltage criterion, c.f. hor- 
izontal dotted lines. The lower graph in Fig. 4(b) shows 
the corresponding Ic{B) measurement for the O-tt junc- 
tion together with a theoretical curve, using 6 ~ 0.33. 
Although the overall agreement between the two curves 
is reasonable, again the shape of the minima is not re- 
produced well. 

To further discuss the observed discrepancies we either 
have to assume, that and are non-uniform over the 
junction length, which would be contradictory to the ob- 
servations at the reference junctions, or we should con- 
sider effects caused by a possible remanent magnetization 
of the F-layer, which can be different in the and tt part, 
plus the possibility that the magnetic fiux generated by 
the applied field may be enhanced by the magnetic mo- 
ment of the F-layer. Also, the effective junction thickness 
A may be different in the and tt parts, causing addi- 
tional asymmetries. To account for these effects, the local 
phases in the two parts may be written as 

fi%x)^cbo + iip%+^°M)x/Lo (3) 
H^{x)^c^o + {^1+^m)x/L^. (4) 

Here, 0o is an initial phase to be fixed when calculat- 



ing the total critical current, (pj^^ are the fluxes, nor- 
malized to $o/27r, that are generated by the (ID) y- 
component of the in-plane remanent magnetizations in 
the and tt parts, respectively, ip^^^ are the normal- 
ized fluxes through the junction generated by the applied 
magnetic field. In the following we parameterize cp^j^'J' as 
= <Pa/(1 ± Sm) and (^^^ as (p^'' = (^s(l ± (5^), re- 
spectively. We further set Lq = L^^ = L/2 which is the 
case for the sample discussed here. 

To obtain the junction critical current as a func- 
tion of the applied magnetic field, wc first calculate the 
currents lojn in the and tt parts via 

.0 

Io= sin(/(x)) dx, 

J-Lo 



It: = jc sin in'^ {x)+TT)dx, 
Jo 

and maximize Jq -I- /,r with respect to (j)Q for each value 
of the applied magnetic field. 

Wc first address the effect of the parameters 6, (pM, 
Sm, and Sb on the I'^~'^{(pb) patterns, c.f. Fig. 5(a) to 
(d). 

If only a jc asymmetry is considered, as shown in Fig. 
5(a), using definitions jj? = jc(l + <^), jc = ic(l — <5), and 
jc = Q.b{jc + jc), one finds that with increasing asym- 
metry i5 the central minimum increases (for 5 ~ I one 
reaches the extremum of a non-stepped junction with 
length Lq, while the tt part becomes "non-Josephson" 
with jc — 0). However in all cases the first side max- 
ima remain symmetric and the side minima reach zero 
current. 

Next we would like to take into account the effect of 
the fiux generated by remanent magnetizations. If we 
consider only a non-zero magnetization, i.e. (pM ^ 0, 
with all other parameters being zero, the Idps) curve 
gets shifted along the field axis, since the total fiux in 
the junction is just the sum of applied field and magne- 
tization. This can be seen in Fig. 5(b) (black curve). By 
adding an asymmetry 6m the side minima get bumped 
and at the same time the maxima decrease (c.f. Fig. 5(b) 
red curve). However the /c(<Ps) curve is still symmetric 
with respect to the central minimum. This changes by 
adding an additional asymmetry (5 7^ in the critical 
current densities. Now the two main maxima get asym- 
metric and the side minima get bumped (blue curve). 

Now we want to consider the effect of asymmetric fiux 
in the and tt halves, i.e. we look at 5b 7^ 0. In Fig. 5(c) 
wc show the results obtained by increasing 6b with the 
other parameters kept at zero. The increase of 6b leads 
to bumped minima and decreased side maxima. The re- 
sulting Ic{<pb) curves looks similar to the ones shown 
in Fig. 5 (b) with asymmetries in the magnetization 6m- 
The comparison reveals that the 6b parameter acts much 
stronger than 6m- The Ic{(Pb) curve is still symmetric 
with respect to the central minimum. 
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FIG. 5; (Color online) Calculated magnetic diffraction pat- 
tern Ic~^{ipB) for a short O-tt junction with asymmetries in 
the critical current densities and in the magnetizations of the 
and TT part: (a) effect of the asymmetry parameter of the 
critical current density 5; (b) resulting pattern with additional 
remanent magnetizations (average value ipM and asymmetry 
(5a/); (c) effect of the asymmetry parameter Sb caused by the 
applied flux; (d) effect of S for nonzero values of </3i\/ — 0.4 
and 5b = 0.06. 
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data 

theory: &,=-1.2 
theory: 5 =-3 




In Fig. 5(d) we add a remanent magnetization without 
asymmetry, i.e. (pM ^ and (5a/ = 0, and allow asym- 
metric critical currents (5 7^ 0. As one can sec the maxima 
remain symmetric whereas the minima get slightly asym- 
metric. 

We further note that the calculated Ic{B) patterns are 
identical if we simultaneously change the sign of S, Sm, 
and 6b- Thus the Ic{B) pattern of the O-tt junction only 
does not allow to identify which parameters belong to the 
and TT part. However the additional information on the 
(temperature dependent) critical current densities of the 
reference junctions may allow a clear identification of 
and TT. 

Using the above findings on the parameters (5, ifAi, 
Sm, and (5b we next discuss our experimental data. For 
the non- vanishing central minimum in Ic{B) a critical 
current asymmetry 6 is required and the shift along the 
magnetic field axis can solely be caused by a finite value 
of ifM- Thus there are only two non-trivial parameters 
(<5m, Sb) left to reproduce the remaining features of the 



FIG. 6: (Color online) Ic{B) patterns of O-tt junction: com- 
parison of experimental data and fitted magnetic diffrac- 
tion pattern /^""(B) for T = 2.65 K using 5 = 0.18 and 
<^A/ = —0.1. In (a) Sm has been varied at fixed Sb = 0, 
in (b) Sb was varied with fixed Sm = 0, and in (c) Sm and Sb 
are varied. 



experimental data. 

If one allows for an asymmetry in the remanent mag- 
netizations only, i.e. Sm 7^ and 6b = 0, it is not possi- 
ble to reproduce the experimental Ic{B) at low and high 
magnetic fields at the same time. The resulting curves 
can be seen in Fig. 6(a). For large 6m = —3 (Fig. 6(a) 
dashed green line) the fit works well for high fields but 
fails in the first side minima. With a smaller value of 
6m = —1.2 (Fig. 6(a) solid red line) the situation is 
opposite. 

By contrast the parameter 6b (with (5j\/ — 0) leads 
to a good agreement between the theoretical and expcr- 
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FIG. 7: (Color online) Ic{B) patterns of O-vr junction: com- 
parison of experimental data and fitted magnetic diffraction 
pattern for T = 4.2K (solid line). The dashed line 

includes the effect of a finite voltage criterion Kr. 



imental Ic{B) pattern. This is shown in Fig. 6(b) where 
wc used 5b — 0.059. There arc only small asymmetries 
near the side maxima and minima that cannot be repro- 
duced for the case 5m =0. If we use both asymmetry 
parameters we get an excellent agreement of the theory 
with the experimental data, as shown in Fig. 6(c) for the 
T = 2.65 K data. 

To further test the fit procedure we now use the T = 
4.2 K data and assume that the magnetic parameters 
remain the same as for T = 2.65 K. By contrast 5 will 
change due to a different temperature dependence of j° 
and j^, as already discussed above. For 5 = 0.33 we get 
a reasonable agreement, as shown in Fig. 7. 

The T = 4.2 K fit is apparently not as good as the 
T = 2.65 K fit. However note that at 4.2 K the detection 
limit is much higher and the minima in Ic{B) are limited 
by the finite voltage criterion. Still some of the bumps 
appear at the same values of applied field both in the 
experimental and theoretical curve. 

For the sake of completeness we also consider the ef- 
fect of th e finite voltage criterion. Using the expression 
V = R\/P — Ic describing the current-voltage charac- 
teristics of a Josephson junction in the framework of 
the resistively shunted junction (RSJ) model ' '' "' we get 

a corrected /^(B) via IcMB) = ^{V^JWTTJbY ^ 
where Ic{B) refers to the theoretical curve (solid curve 
in Fig. 7) . The corrected curve is shown as dashed (green) 
line in Fig. 7. As can be seen the data are reproduced 
perfectly. 



Discussion of the obtained parameters 

Finally, we would like to discuss the parameters which 
arc obtained by fitting the experimental data. As already 
mentioned above, the parameters 5, (pM, 5m, and Sb al- 



low to find the different parameter sets for the two halves 
of the junction. For the distinction between "0" and "tt" 
additional information is needed, which we get from the 
reference junctions. 

The parameter 5 allows to extract the absolute val- 
ues of the critical current densities in the two parts. We 
get an almost temperature independent critical current 
density in the one half with ji(4.2 K) « jc^(2.65 K) « 
2.3 A/cm^. By contrast, for the other part we find a tem- 
perature dependent current density with jc(4-2 K) « 
1.2 A/cm^ to j^(2.65 K) w 1.6 A/cm^. A comparison 
with the temperature dependencies of the reference junc- 
tions allows the identification that the first part has to be 
coupled, whereas the second part is tt coupled. The ab- 
solute values of jj? of the and O-tt junction are approxi- 
mately the same whereas jj! of the O-tt junction is reduced 
by w 0.5 A/cm^ as compared to the tt reference junction, 
although the temperature dependence looks very similar. 
This indicates a slightly reduced thickness of the F-layer 
in the tt part of the O-tt junction, c.f. Figure 1(b). Taking 
the data of Ref. 8 the difference in thickness can be esti- 
mated to be « 0.7 A. This may be caused by some gra- 
dient of the ferromagnetic thickness along x direction on 
the chip, as the distance between reference and stepped 
junctions on the chip is about 2 mm. 

The parameters related to a different remanent mag- 
netization in the and tt part, i.e. (pM = —0.1 and 
5m = —0.3, seem reasonable. The magnetization is of 
the order of 10^"^ of a fully saturated magnetization, in- 
dicating that the F-layer is in a multi-domain state. Note 
that the resulting magnetization of the tt part is larger 
than the magnetization of the part, which seems real- 
istic due to a thicker F-layer in the tt part. In fact, the 
ratio of the F-layer thicknesses d2/di is very close to 1, 
so, assuming that magnetization is proportional to the 
volume of the F-layer in each part, it is quite difficult 
to explain the above value of However, if one as- 
sumes that there is a dead layer of thickness ddead one 
can calculate its value from 

dl — ddoad _ 1+ 5m 
d2 — rfdcad 1 — 5m 

to be ddcad ~ 4.7 nm. This value is somewhat larger 
than ddcad ~ 3.1 nm estimated earlier from a jcidp) fit 
made for a different run of the same fabrication process'"'. 
However, as we see from Figs. 6 (b) and (c) the change 
in 5m from to —0.3 affects only the tiny features on 
the Ic{B) curve. Thus, the value of 5m cannot be found 
from this fit very exactly. 

Besides the current asymmetry 5 the most important 
parameter for our experiment is the asymmetry param- 
eter 5b- Using a finite 5 and 5b = 0.059 almost all fea- 
tures could be reproduced very well. The addition of 
the parameters related to remanent magnetizations lead 
to minor improvements in the agreement of theory and 
experiment. In the following we want to discuss three 
possible scenarios causing the asymmetry 5b- 

First, the effect could be caused simply via the fab- 
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rication procedure of the junction. In the part of the 
junction the SF bilayer was deposited in situ whereas the 
Nb cap layer in the tt part was deposited after an etch- 
ing process. Thus the properties, such as the mean free 
path and hence the London penetration depth A, of the 
Nb cap layers in the two halves could easily differ by few 
percents. 

Second, one could think of a paramagnetic component 
in the magnetization. As already discussed above, the 
F-layer is expected to be in a multi-domain state with 
a small net magnetization in-plane. An external field 
applied in-plane could cause a reconfiguration of the do- 
mains. In the two halves the pinning of the domains may 
be different, due to the different thicknesses and the dif- 
ferent treatment. This would result in a asymmetric field 
dependent magnetization. 

A third possibility is the appearance of an enhanced 
flux penetration due to inverse proximity effect, causing 
a correction in the London penetration depth. Due to 
the reduction of the order parameter in the vicinity of 
the ferromagnetic layer, the effective penetration depth 
might be enlarged. In order to estimate this effect, we cal- 
culated numerically the space-dependent superfluid den- 
sity ns{z) in the superconducting and the ferromagnetic 
part of the SF bilayer using the quasiclassical approach" '. 
Herein we used the parameters of our SF bilayer, which 
were already obtained in Ref. 24 by fitting the experi- 
mental data of Ref. 8 . By using the (London) expression 
A(z) oc ns{z)~^-^ we obtained the spatial dependence of 
the penetration depth. Then we used the second London 
equation V'^B{z) = B{z)/}?{z) to calculate the magnetic 
field B{z) numerically. We define the effective penetra- 
tion depth as \L,eS = Ai^es/^, with 4>eff and $ being 
the flux in our SF bilayer with and without inverse prox- 
imity corrections. For our SIFS junctions with a thick- 
ness dp Ri 5 nm of the ferromagnet and t2 = 400 nm of 
the top electrode we get Xl,cS = I.OOSAl at T = 2.65 K. 
Therefore in our case the inverse proximity corrections 
are negligible. In addition the corrections due to inverse 
proximity effect would be opposite in sign, i.e. Sb < 0, 
in contrast to Sb = -f 0.059 found for our junction. 

By looking at the other two scenarios it seems natural 
that the fabrication procedure causes the observed Sb 
asymmetry. However at the moment, we cannot exclude 
a field-dependence of the magnetization. A clarification 
deserves further investigations. 

IV. CONCLUSIONS 

In this paper we presented a detailed analysis of the 
magnetic field dependence in the critical current, Ic{B), 



in 0, TT, and O-tt SIFS Josephson junctions. The length of 
the junctions is smaller than the Josephson length. The 
Ic{B) pattern of the and the tt junction can be well 
described by the standard Fraunhofer pattern, valid for 
a homogenous, short junction. The central maximum of 
this pattern is typically shifted from zero by some percent 
of one flux quantum, pointing to a weak in-planc magne- 
tization of the F-layer. The magnetization is of order of 
10^'^ of a fully saturated magnetization, indicating that 
the F-layer is in a multi-domain state. 

The Ic{B) pattern of the O-tt junction exhibits the cen- 
tral minimum, well known for this type of junction. How- 
ever the critical current at this minimum is non-zero, 
pointing to an asymmetry in the critical current densi- 
ties in the two halves of the junction. In addition Ic{B) 
exhibits asymmetric maxima and bumped minima that 
cannot be described exclusively by critical current asym- 
metries. A detailed explanation of these features requires 
the consideration of asymmetric fluxes generated in the 
and TT parts of the junction. A careful analysis of the ex- 
perimental data and our model showed that the majority 
of the observed discrepancies are due to a field-dependent 
asymmetry of the fluxes in the and tt part. The ef- 
fect could cither be caused by a small, fleld-dependent, 
in-plane magnetization of the F-laycr or by a difference 
in the penetration lengths, which most naturally can be 
due to the fabrication technique. In principle, this ef- 
fect should also be present in the /c(i3)'s of the reference 
junctions. However, here the effect only leads to a small 
scaling factor for the magnetic field, which is too small 
to be detectable in experiment, e.g. if the effects of field 
focusing are considered. 

The model discussed in this paper on the basis of O-tt 
junctions can be extended, e.g., to SIFS junctions hav- 
ing step-like jc{x) profile' '. or laterally ordered ferromag- 
netic domains. 
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